Distinct Pairing Symmetries in Ndi^Ceo^CuO^-y and Lai.sgSro.nCuO^ Single Crystals: 
Evidence from Comparative Tunnelling Measurements 
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We used point-contact tunnelling spectroscopy to study the superconducting pairing symmetry of 
electron-doped Ndi.ssCeo.isCuOi-y (NCCO) and hole-doped La 1 .sgSro.uCu04 (LSCO). Nearly 
identical spectra without zero bias conductance peak (ZBCP) were obtained on the (HO) and 
(100) oriented surfaces (the so-called nodal and anti-nodal directions) of NCCO. In contrast, LSCO 
showed a remarkable ZBCP in the nodal direction as expected from a d-wave superconductor. 
Detailed analysis reveals an s-wave component in the pairing symmetry of the NCCO sample with 
A/ksTc = 1.66, a value remarkable close to that of a weakly coupled BCS superconductor. We 
argue that this s-wave component is formed at the Fermi surface pockets centered at (±7r,0) and 
(0,±7r) although a d-wave component may also exist. 

PACS numbers: 74.50.+r, 74.72. Jt, 74.45. +c 



I. INTRODUCTION 

For hole-doped cuprate superconductors, there is a 
well-known electronic phase diagram characterized by a 
dome-like superconducting region when p is above a cer- 
tain threshold It was found that @ instead of dop- 
ing holes into the Cu-0 plane, one can also achieve su- 
perconductivity by doping electrons into the Cu-0 plane 
in systems like Ln2_ x Ce2;Cu04 with Ln = Nd, Pr, La, 
etc. One thus naturally questions whether the supercon- 
ducting mechanism is the same between the hole- and 
electron- doped cuprates. Among many superconduct- 
ing properties, the symmetry of the order parameter is 
an important one. While the symmetry of the supercon- 
ducting order parameter is believed to be d-wave for the 
hole-doped region the situation for electron-doped 
materials is highly controversial. For example, Angle- 
resolved photoemission spectroscopy (ARPES) [JL spe- 
cific heat [5], phase-sensitive scanning SQUID 6], bi- 
crystal grain-boundary Josephsoniunction[7| and some 
penetration depth measurements [g, Q indicate a d-wave 
symmetr y. I n addition, Raman scattering [To| and recent 
ARPES |ll| experiments show a nonmonotonic d-wave 
order parameter. However, this has been contrasted by 
tunnelling [l2|, EH and some other specific heat [l4[ and 
penetration depth measurements [3, EH, [TtJ • In partic- 
ular, there may be a crossover from d-wave to s-wave 
symmetries by changing the doped electron concentra- 
tion [H, EH or decreasing temperature [20| . Although 
such a crossover may explain the conflicting experimen- 
tal results on the pairing symmetry in the electron-doped 
cuprates, the characteristics of such possible s-wave pair- 
ing symmetry has yet to be established. 

In this paper, we report directional tunnelling mea- 
surements on single crystals of optimally electron-doped 
cuprate Ndi^Ce^.i^CuO^-y (NCCO). By injecting cur- 
rent along either the Cu-Cu bond (110) or Cu-0 bond 



(100) direction, we obtain nearly identical tunnelling 
spectra indicating an s-wave component of the pairing 
symmetry in this material. For comparison, similar mea- 
surements were carried out on underdoped p-type single 
crystals of Lai. ^StquCuOa (LSCO), and we observe 
clear zero bias conductance peaks (ZBCP) in tunneling 
spectra along the (110) direction as expected from the 
d a; 2_ 1 ,2-wave symmetry. Our results thus indicate that 
the optimally doped NCCO has at least an unavoidable 
s-wave component, which is in contrast with the case in 
hole-doped LSCO where a pure d-wave has been estab- 
lished. 



II. EXPERIMENTAL DETAILS 

We grew the NCCO and LSCO single crystals using 
the travelling-solvent floating-zone technique [U • As 
shown in Fig. [T] the resistive curve measured on the 
NCCO sample indicates a zero-resistance temperature 
at 25. IK, AC susceptibility shows the onset of bulk su- 
perconductivity at T c w 25. 6K. The LSCO sample has 
a T c w 28K characterized by AC susceptibility mea- 
surement. The single crystals used in the tunnelling 
measurements were cut into rectangular flakes of size 
2.5 x 2.5 x 1mm 3 for NCCO and 1.5 x 1.5 x 0.8mm 3 
for LSCO with the long axis along a (100) or b (010) and 
short axis along c (001), one corner of the crystal was 
then cut off to expose (llO)-oriented crystal planes. The 
directional point-contact tunneling measurements were 
carried out by pointing a Pt/Ir alloy tip or Au tip to- 
wards the specified directions as shown in the left inset 
of Fig. [T] The tip's preparation and the equipment de- 
tails are described in Ref.[23[. In order to reduce quasi- 
particle scattering in the barrier layer and hence obtain 
high quality data, the nonaqueous chemical etch was used 
to attenuate the insulating layer on the sample surface 
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FIG. 1: The temperature dependence of resistivity for the 
sample NCCO. The insets show the AC susceptibility of 
NCCO and LSCO. The schematic diagram of the point con- 
tact configuration are presented in the upper left of the main 
panel. 



immediately before mounting the sample on the point 
contact device. Since it takes about 20-40 minutes to 
insert the sample mounted device into the Helium gas 
environment in the dewar, the sample surfaces were ex- 
posed to air during this period. For chemical etching, we 
use 10% HC1 in absolute ethanol for several seconds to 
several tens of seconds for NCCO and 1% Br 2 in absolute 
ethanol for several minutes to ten minutes in the case of 
LSCO. Typical four-terminal and lock- in technique were 
used to measure the I ~ V curves and the differential 
resistance dV/dl vs V of the point contacts. 



III. RESULTS AND DISCUSSIONS 

A. Tunneling spectra of Ndi.gsCeo.isCuOi-y 
(NCCO) 

Fig. [2] shows the raw data of the conductance (a) of 
Au/NCCO point contact for various temperatures from 
2K to 30K. In order to show the electronic stability of the 
experiments, we simultaneously present in Fig. [3](a) the 
a ~ V curves obtained from the lock-in technique and 
DC I ~ V measurement, respectively. It is obvious that 
these two curves merge into each other although the lat- 
ter one has a higher noise. We also compared the a ~ V 
curves recorded by both positively and negatively scan- 
ning the bias voltage. As shown in Fig.^b), the complete 
overlap of the data with different scanning directions in- 
dicates the good thermal stability during measurements. 

It is well known that the high T c superconducting 
cuprate compounds react readily with atmospheric H2O 
and CO2 to form insulating hydroxides and carbonates 
on the surface and at grain boundaries [24j |. While such 
poor conductive layer may be a natural barrier needed 
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FIG. 2: Raw data of the directional spectral measurements, 
(a) The temperature dependent tunneling spectra measured 
along (100) direction. The curves have been shifted for clarity. 
The temperature increases from the bottom up with the steps 
of IK (from 2K to 22K) and 2K (from 22K to 30K). The thick 
solid line denotes the data at 26K which is around T c . (b) An 
illustration of constructing normal conductance background 
according to the functional form of the spectra above T c . (c) 
Normalized spectrum for (110) direction at T = 2K. 
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FIG. 3: Comparison between the a ~ V curves (a) obtained 
from lock-in technique (solid line) and DC I ~ V measure- 
ments (open circles); (b) measured with positive (solid line) 
and negative (open circles) bias scanning. 
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FIG. 4: The normalized a ~ V curves obtained by differ- 
ent surface treatments and using different tips. Note that 
the spectra measured on the post-etched surfaces are much 
sharper than that of pre-etched ones, accompanied by a re- 
markable decrease in the junction resistance from several hun- 
dreds of Ohms to below 10 Q as shown in Table [I] 



TABLE I: List of the experimental conditions for the different 
spectra presented in Fig. [4] 



Label 


Chemical Etch 


Tip 


R20mv (n) 


(110)-a 


None (As cut) 


Pt/Ir Alloy 


197.0 


(110)-b 


By 10% HC1 (20s) 


Pt/Ir Alloy 


75.7 


(110)-c 


By 10% HC1 (60s) 


Au 


5.3 


(100)-a 


None (As cut) 


Pt/Ir Alloy 


170.0 


(100)-b 


By 10% HC1 (20s) 


Pt/Ir Alloy 


19.9 


(100)-c 


By 10% HC1 (60s) 


Au 


2.9 



in tunneling experiments, they may introduce large scat- 
tering factor (just like the situation of fabricated planar- 
junction) and thus result to less sharpening spectra than 
those measured by the vacuum barrier based STM [2j|. 
In order to reduce the scattering factor and obtain sharp 
spectra, we used no naq ueous chemical etch to treat the 
sample surfaces [HI, HH]. In Fig. [4j we present the mea- 
sured spectra on the surfaces as cut, after the first etch 
and after the second etch, all the spectra have been nor- 
malized for comparison. The detailed experimental con- 
ditions are listed in Table HI The experimental results 
indicate that the chemical etch can effectively reduce the 
thickness of the barrier layer on the surface of NCCO 
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FIG. 5: The normalized a ~ V curves measured at a fixed po- 
sition on the sample surface with various junction resistance. 



single crystal and hence decrease the scattering factor. 

We have also studied the dependence of the measured 
spectrum on the junction resistance which can be con- 
trolled by adjusting the tip-to-samplc distance with a 
differential screw. Fig. O shows the spectra measured at 
a fixed position on the sample surface as a function of 
junction resistances, in these measurements the sample 
surfaces had been chemically etched. When the junc- 
tion resistance decreases from several tens of Ohms to 
several Ohms, the spectral shape becomes sharper (refer 
to Fig. Ufa)). For small resistance, the spectral shape 
changes hardly with the resistive values and no Andreev 
reflection-like spectrum occurs [27[ (Fig. [5]) . This indi- 
cates that although the surface contaminant phase has 
been reduced, it still exists as a solid barrier layer due 
to the exposure to air. After the metal tip reaches the 
sample surface, the barrier layer is abraded at first and 
its thickness decreases with the increasing pressure. Con- 
sequently, the measured spectrum becomes sharper due 
to the weakening of the inelastic scattering of the inject- 
ing quasiparticles near the normal-metal/superconductor 
(N/S) micro-constriction. Further pressure of the tip on 
the sample surface may simply flatten the point, giving 
a larger contact area over the same minimal thickness of 
a tenacious barrier layer which can not be entirely flaked 
off |28j . In this case, the junction resistance becomes 
still smaller while with no obvious change on the mea- 
sured spectrum until the junction is damaged eventually. 
Therefore, the none-zero zero bias conductance should 
be ascribed to the stronger scattering effect on injecting 
quasiparticles and lower height of the surface barrier than 
that of the vacuum barrier. Nonetheless, such surface 
contaminants resulting from air exposure form orders of 
magnitude slower on NCCO than on many hole-doped 
cuprates which results from the absence of reactive alka- 
line earth elements in NCCO [26| . This is also the reason 
why the measured spectra of LSCO are poorer than that 
of NCCO using the same method, as shown in the section 
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FIG. 6: The magnetic field dependence of the tunneling spec- 
tra measured at (a) 4K and (b) 10K. The very small modifi- 
cation of the spectral shape in field up to 7T indicates that 
we indeed measured the bulk superconductivity of NCCO. 



III-C in this paper. In the following, we will focus our 
attention on the spectra with the resistance from 1 to 15 
f2, namely, in the regime with good spectral stability and 
weaker scattering effect. 

In MaglabExa-12, the magnetic field can be applied 
along the tip's direction, namely, parallel to the ab- 
plane. The field dependence of the tunneling spectra 
were measured at various temperatures from 4K to 16K. 
The changing tendency of the spectral shape with in- 
creasing field is identical at all measured temperatures. 
The typical results of 4K and 10K in Fig. [5] suggest that 
the field induced smearing of the spectrum is very slight, 
which is consistent with the much higher H C 2 than 7T 
for H || ab — plane. 

To sum up the main characteristics of the measured 
spectra, we note that all the spectra are near the tun- 
neling limit with two clear coherence peaks at symmetric 
positions of the bias voltage (much sharper than that of 
the previous works studied by both STM [l3|, [29[ and 
point contact [30j|) with no Andreev reflection- like spec- 
trum. Second, the spectral shape is nearly identical for 
the (110) and (100) directions. Thirdly, there was no ev- 
idence for zero-bias conductance peak (ZBCP) expected 
in the (110) (or the so-called nodal direction) of the d- 
wave superconductors [U, H3, Hj|. Here we emphasize 
that the results are independent of the positions on the 
sample surface. 
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FIG. 7: The comparison of fitting the normalized conductance 
measured along (llO)-direction at 2K to the s-wave model and 
d a! 2_j / 2-wave model. For, d-wave fitting, a denotes the angle 
between the normal direction of the sample surface and the 
crystallographic axis a. Note that the s-wave fitting is much 
better than the d-wave fitting. 



B. Theoretical model - s-wave BTK theory 

As discussed above, if the superconducting NCCO 
has d-w&ve pairing symmetry, a ZBCP should be ob- 
served in the (llO)-oriented tunneling spectrum. In other 
words, our experimental results suggest that the opti- 
mally doped NCCO has a s-wave symmetry. In or- 
der to explain this viewpoint, we present in Fig. [7] the 
best fitting of the (llO)-oriented tunneling spectrum with 
the formulas corresponding s-wave and cZ-wave models 
respectively. For the simulations, the extended BTK 
model HEHl! was accepted by selecting a constant gap 
value for s-wave symmetry and the anisotropic gap of 
A(#) = Aqcos(29) for d x i_ y i symmetry, where 9 is the 
polar angle measured from the crystallographic axis a. In 
this model, two parameters are introduced to describe the 
necessary physical quantities, i.e., the effective potential 
barrier (Z) and the superconducting energy gap (A). As 
an extension, the quasiparticle energy E is replaced by 
E — iT, where T is the broadening parameter characteriz- 
ing the finite lifetime of the quasiparticles due to inelastic 
scattering near the N/S micro-constriction [34], l3q . In 
a real N-I-S junction configuration, total tunneling con- 
ductance spectrum includes the integration over the solid 
angle. In the case of two dimension, it reduces to the in- 
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FIG. 8: (a) The normalized spectra measured at a fixed po- 
sition on the sample surface with various junction resistance 
in a large resistive range. The open circles denote the ex- 
perimental data and the solid lines indicate the fitting to the 
s-wave extended BTK model, (b) The fitting parameters for 
the spectra presented in (a). 



tegration over the injection angle from — ir/2 to tt/2, as 
done in this work. Fig.[7]Jb) clearly shows that the <i-wave 
theoretical simulation along the (110) direction deviates 
from the experimental data. When the normal direction 
of the sample surface departs a small angle from the crys- 
tallographic axis a, a ZBCP will appear, accompanied 
by the remarkable depression of the coherence peaks. As 
shown in Fig. EJa) , even the d-wave simulation absolutely 
along (100) direction can not fit the experimental data in 
the range below superconducting gap. However, the cal- 
culated curves in terms of the s-wave theory are in good 
agreement with the experimental results both in (100) 
and (110) directions. 

It should be pointed out that the values of Z and T 
are related to the selected form of the background (such 
as the slope coefficient of the high bias straight-line seg- 
ments) in a certain extent, which may be the main ar- 
tificial uncertainty in the data analysis. So we did not 
focus our attention on the absolute values of these pa- 
rameters or the quantitative comparison between the dif- 
ferent spots. However, it is also found that the energy 
gap A is insensitive to the selection of different back- 
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FIG. 9: The temperature dependent normalized conductance 
and the s-wave BTK simulations for (a) (110) direction and 
(b) (100) direction. All the curves have been shifted upwards 
except the lowest one for clarity. The temperature increases 
from 2K to 20K with a interval of IK (from the bottom up). 
The theoretical simulations are denoted by red solid lines. 



ground in the fitting process. Moreover, when we chose 
a general procedure to construct the background for all 
spectra measured on a fixed spot on the sample surface, 
the variation of the fitting parameters of Z and T should 
be physically meaningful. As an example, we presented 
in Fig. Efa) the normalized spectra measured at a fixed 
spot with different junction resistance, which have been 
theoretically fitted to the s-wave BTK model and the fit- 
ting parameters are shown in Fig.[8jb). It is noted that 
the derived energy gap is nearly constant for all spectra, 
while the barrier height Z and broadening parameter T 
continuously decreases with the decreasing junction resis- 
tance. This is physically reasonable and in good agree- 
ment with the foregoing discussions, namely, the abrasion 
of the contaminant layer will depress the barrier height 
accompanied by the weakening of the quasiparticle scat- 
tering near the N/S micro-constriction. 

The good spatial repeatability of the tunneling spectra 
allowed further investigations on their temperature de- 
pendence. The spectral shape change continuously until 
above 24K (T c ~ 26K), thus suggesting that the tun- 
neling spectra reflect the bulk superconductivity of the 
NCCO sample (Fig. [2Ja)). The temperature dependent 
normalized spectra for both (110) and (100) directions 
are presented in Fig. accompanied by the s-wave BTK 
fitting curves denoted by the red solid lines. Because 
the spectra were affected by the critical current effect at 
higher temperature near T c 36], it is difficult to deter- 
mine the background for these temperatures. Therefore, 
these results have not been included in our theoretical 
analysis. All the parameters of the simulations shown in 
Fig. [9] are presented in Fig. [10] (the barrier height Z and 
the broadening parameter T) and Fig. [TT] (superconduct- 
ing energy gap A). The junction resistances measured at 
20mV (Riomv) are also given in Fig. [TO] 
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FIG. 10: The junction resistance measured at 20mV, the best 
fitting parameters of barrier height Z and broadening param- 
eter r corresponding to the simulations shown in Fig. O The 
junction resistances RwmV of different spots on the sample 
surface possess the similar decay law with increasing tem- 
perature, indicating that the slight thermal expansion of the 
point-contact device. Correspondingly, the values of Z and V 
decrease with increasing temperature at first, then tend to be 
constant when the junction resistance becomes small enough 
because the minimal thickness of the barrier layer is achieved 
and the point flattening occurs. 



1.0 



g 0.8 



< 0.6 



0.4 



NCCO 



0.0 




Direction A ( meV 
□ (100) 3.55 
O (110) 3.67 
A (110) 3.65 
BSC theory 



0.2 



T/T„ 



FIG. 11: The determined A(T) relations by fitting the tem- 
perature dependent tunneling spectra to the s-wave BTK the- 
ory. The solid line indicates the BCS prediction. 



The magnitude of A = T/A(0) < 20% presented is 
larger than that of the clean point-contact junction be- 
tween normal metal and conventional superconductors 
which in general is smaller than 10%. However, con- 
sidering the much shorter coherence length and faster 
exterior degradation of cuprates than that of the conven- 
tional superconductors, one can easily understand the 
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FIG. 12: The point contact spectrum measured at 2K on 
the Nb foil exposed in air for long time. The solid line 
denotes the s-wave BTK fitting. The obtained value of 
A — r/A(0) = 0.28 is much larger than that obtained from 
Au/NCCO tunneling in this work. 



larger A. Actually, a large A (> 20%) is often obtained 
on the oxidized surface of Nb foil, a typical conventional 
superconductor, as shown in Fig. [T3J The lowest value 
of A = 17% achieved in this work is much smaller than 
the previous report of A 1 by STM measurements on 
NCCO single crystal 0. 

As shown in Fig. llH the normalized gap function A(T) 
determined by the fitting procedure can be well described 
by BCS theory, yielding a gap ratio A(0)/fcsT c « 1.66, 
which is very close to the theoretical value of 1.76. Such 
consistency between the tunneling data of cuprates and 
the theoretical model over a wide temperature range (be- 
tween OK and T c ) is surprising. In order to determine 
the repeatability of such temperature dependent mea- 
surements, we plot in Fig. [11] two A(T) relations ob- 
tained at two different positions for the (110) direction. 
All the A(T) functions obtained along both (110) and 
(100) directions and from different positions follow the 
BCS theory in a normalized scale. The T c values derived 
from this figure is between 24.5K~26K, near the bulk 
transition temperature T c w 2b.6K. 

The distribution of the superconducting energy gap 
A(0) derived from the spectral measurements at many 
different positions are shown in Fig. 1131 There is no obvi- 
ous difference between the data obtained from (110) and 
(100) directions. The value of A(0) varies in a narrow 
range from 3.50 to 3.70, indicating the good homogene- 
ity of the superconductivity in the investigated regions. 



C. Tunneling spectra of Lai.ggSro.nCu04 (LSCO) 

As a comparison, the directional tunneling spectra 
along both (100) and (110) directions were also studied 
on the LSCO single crystal, which has a typical d x 2_ y 2- 
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FIG. 13: Statistical chart of the A(0)-values determined by 
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surfaces. The y-axis denotes the occurring times for a given 
gap- value specified by the ir-axis. 
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FIG. 15: Normalized tunneling spectra of T = 2K. Thick and 
thin solid lines indicate (110) and (100) directions, respec- 
tively. The remarkable difference between the cases along the 
two axis is in good agreement with the ci-wave theory. 
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FIG. 14: Temperature dependent tunneling spectra of 
Lai.89Sro.nCu04 single crystal measured along (a) (100) di- 
rection and (b) (110) direction. The distinct spectral charac- 
teristics of these two directions are obviously different from 
the case of NCCO. 



wave pairing symmetry (as checked recently by specific 
heat [371]). The measured spectra at various tempera- 
tures are presented in Fig. Q3J For the (100) direction 
at low temperatures, two clear coherence peaks appear 
in the spectra and no ZBCP can be observed. While for 
the (llO)-oriented spectra, a prominent ZBCP was ob- 
served, accompanied by the disappearance of the coher- 
ence peaks. With the increasing temperature, all these 
characteristics become weaker and weaker and eventu- 
ally vanish around T c . As mentioned above, photoemis- 
sion experiments have proved that surface degradation 
resulting from air exposure form orders of magnitude 
faster on LSCO than NCCO, due to the existence of 
reactive alkaline earth elements in LSCO [26|. There- 
fore, the surface of LSCO degraded quickly before being 
put into the helium gas environment, which directly af- 
fected the subsequent spectral measurements. That is, 
the injecting quasiparticles felt strong scattering effect 
from the surface barrier layer (resulting in a very large 
spectral broadening factor), i.e., the measured spectra 
were badly smeared. As shown in Fig. 1141 the super- 
conducting characteristics nearly completely disappear 
above 20K which is much lower than the bulk critical 
temperature T c w 28K. 

In order to further demonstrate the spectral differences 
along the two directions, the spectra of T = 2K have 
been normalized using the higher bias backgrounds, as 
shown in Fig. [151 Although the strong scattering effect 
has severely smeared the spectra and made the quanti- 
tative analysis difficult, the coherence peaks around gap 
energy are depressed and a ZBCP appears when the tun- 
neling direction changes from (100) to (110) as expected 
from the d-wave theory. These spectra are also consis- 
tent with the results obtained on the La2_ a; Sr z ,Cu04/Ag 
junctions fabricated using a ramp-edge technique [381 ]. 
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TABLE II: List of the parameters for the fitting presented in 
Fig. m 



Label 


A (meV) 


Ai (meV) 


Z 


T (meV) 


r/A 


(a)-(110) 


3.58 


0.46 


2.95 


0.68 


0.17 


(a)-(100) 


3.45 


0.44 


2.80 


0.84 


0.21 


(b)-(HO) 


4.50 


None 


3.20 


0.54 


0.12 


(b)-(100) 


4.30 


None 


3.10 


0.64 


0.15 



D. Theoretical model - Anisotropic s-wave and 
two-band model 

In this section, we try to explain our experimental 
data with other possible models. In Ref . [13j , the au- 
thors ascribed the unphysically large A (A — T/A ^ 1) 
to the unreasonable assumption of isotropic s-wave sym- 
metry. They found that if the maximum value of T/A 
was assumed to be 0.2 similar to YBCO, reasonable fit- 
ting is obtained by assuming anisotropic s-wave sym- 
metry, namely, A(0) — Ao + Aicos(40). Comparing 
this previous report with present work, we find that the 
large value of A is mainly due to the surface degra- 
dation. Nonetheless, the anisotropic s-wave symmetry 
seems to be reasonable considering the crystallographic 
symmetry and the topology of Fermi surface (FS) of 
NCCO [Hi, HJ. In Fig. [H|a), the best fitting to such 
anisotropic s-wave model is presented for both (110) and 
(100) directions. Recently, the ARPES experiment re- 
vealed another form of the gap function described as 
A(0) = A o -|cos(20)-O.3cos(60)|. Such A (0) function has 
also been tried in our fittings as shown in Fig. [LoTb) . All 
the parameters for the theoretical simulations of Fig. [T|)] 
are listed in Table |TTJ It is found that the isotropic s- 
wave is the best candidate among the three s-wave mod- 
els mentioned above. 

Recent ARPES experiments [4l|, [42j revealed the 
doping evolution of the Fermi surface (FS) in 
Nd2- a; Ce a ;Cu04. At low doping, a small Fermi pocket 
appears around (7T,0). Upon increasing doping, another 
pocket begins to form around (7r/2,7r/2) and eventually 
at optimal doping x = 0.15 several FS pieces evolve into a 
large curve around (7r,7r). These findings were effectively 
described as a two-band system (ioj . Most recently, Luo 
et al. used a weakly coupled two-band BCS-like model 
to account for the low energy electromagnetic response of 
superconducting quasiparticles in electron-doped materi- 
als [43[ . The special angle dependence of the energy gap 
observed in Ref.[llT| (as shown in the inset of Fig. [T5f b)) 
may be a reflection of the two-band model. In the anal- 
ysis in Ref. [43], both the pairing symmetries of band-1 
and band-2 are assumed to be the d x 2_ y i type in or- 
der to achieve the best fitting for the superfluid density 
data, in which the labels 1 and 2 represent the bands 
contributing to the FS centered at (±7r,0), (0,±7r) and 
(±7r/2,±7r/2), respectively. However, as pointed out by 
the authors, there is a finite excitation gap in band-1 




1m ^10 10 20 

Bias voltage ( mV ) 



FIG. 16: Fitting the directional tunneling spectra to the 
anisotropic s-wave model with gap functions of (a) A(9) = 
Ao + Aicos(46>) and (b) A(0) = A • \cos(26) - O.3cos(60)|. 
The insets are the schematic drawing of the angle dependent 
energy gaps. The best fitting to the first type of anisotropic 
s-wave model has the value of Ai/Aq smaller than 15% which 
is negligible within the fitting precision. If the second model 
is accepted, the simulations are also not bad as a whole, but 
can not fit very well the spectral shape around zero bias. 



since the nodal lines do not intersect with the FS of that 
band if the system is not heavily overdoped. This indi- 
cates that the superconducting state of NCCO is actu- 
ally a mixture of d-wave and s-wave-like pairing states. 
In other words, this analysis can not definitely determine 
whether the pairing symmetry of band-1 is s-wave or d- 
wave type. Based on this two-band model, we calculated 
the directional tunneling spectra by assuming different 
pairing symmetries (s-wave or d-wave) for band-1 while 
definite d x 2_ y 2 for band-2. In these calculations, the ratio 
of the contribution from band-1 to that from band-2 was 
chosen according to the discussions in Ref. [HI]. If only 
the d-wave symmetry of band-1 is accepted, a prominent 
ZBCP will appear in the spectra along (110) direction. 
That is, the s-wave symmetry of band-1 is required in 
the fitting procedure. Moreover, the best fitting needs 
a much larger value of A = T/A for band-2 than that 
for band-1 as shown in Fig. [T71 possibly due in part 
to the strong depression of the d-wave superconductiv- 
ity on the sample surface while the essential origin is yet 
to be found. Such mixture of superconductivity com- 
ing from two bands may be another possible reason of 
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^20 -10 10 20 

Bias voltage ( mV ) 



FIG. 17: Fitting the directional tunneling spectra to the two- 
band model, (a) Along (100) direction; (b) Along (110) direc- 
tion. All the fitting parameters are listed in the correspond- 
ing panels. The angle dependent gap amplitude and sign are 
described in the inset. Dotted lines indicate the s-wave com- 
ponent while the solid lines indicate the d-wave component. 

the contradicting reports on the pairing symmetry by 
different experiments. For example the phase sensitive 
measurements @, 0] may selectively detect the gap in- 
formation from band-2 which crosses the Fermi surface 
near (7r/2,7r/2). In addition, such inter-band mixture 



may also be responsible to the doping dependent pairing 
symmetry observed in Pri- x Ce x CuO^ thin films [19j . 
In any case, the s-wave pairing symmetry appears to 
be an important component in the superconductivity of 
optimally-doped NCCO. 

In summary, the directional tunnelling spectra along 
(110) and (100) axis on the NCCO and LSCO single 
crystals illustrate clearly distinct pairing symmetries. In 
contrast to the results of LSCO (x=0.11), no ZBCP was 
observed for NCCO (x=0.15) along the two different di- 
rections while sharp coherence peaks are existent for both 
directions, which disagrees with the pure d-wave pair- 
ing symmetry. The almost identical spectral shapes for 
the two directions on NCCO can be understood in the 
framework of s-wave BTK theory, leading to a BCS type 
temperature dependence of energy gap with the ratio of 
A/kBT c ss 1.66. The present work provides evidence for 
the s-wave component in the superconductivity of the 
optimally-doped NCCO, which should mainly come from 
the band crossing the FS around (±tt,0) and (0,±7r) if a 
two-band model is accepted. 
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